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1. Introduction 

In the past five years a lot of progress has been achieved in calculating QCD Equation of State 
(for recent reviews see ^, In the most recent calculations the Equation of State (EoS)has 
been evaluated for 2+1 flavor QCD, /. e. in QCD with one strange quark and two light {u,d) quarks 
using various improved staggered fermion actions ^, |6[ |7|]. The most extensive calculations of 
the EoS have been performed with p4 and asqtad staggered fermion formulations on lattices with 
temporal extent N-^ = 4, 6 [§, ^ and 8 [0]. These actions improve both the taste symmetry of the 
staggered fermions as well as the quark dispersion relations. The latter insures that thermodynamic 
observables are /^{a^) improved at high temperatures and thus have only a small cut-off depen- 
dence in this regime. The stout-link action, which has been used for the calculation of the EoS 
on lattices with temporal extent Nr = 4, 6 [Q], only improves the taste symmetry of the staggered 
fermions and therefore has the same large discretization errors at high temperatures as the standard 
staggered fermion formulation. On the other hand the stout action has better taste symmetry which 
is important at low temperatures. 

While at high temperatures the masses of the relevant degrees of freedom, quarks and gluons, 
are small compared to the temperature scale, this is not the case at low temperatures and in the tran- 
sition region. One thus may expect that at these temperatures thermodynamic observables are more 
sensitive to the quark masses, which control the mass of the light pseudo-scalars and eventually are 
responsible for the occurrence of a true phase transition in the chiral limit. Calculations with p4 
and asqtad actions have so-far been performed using light quark masses (m/) which are one tenth 
of the strange quark mass (m^) and correspond to a lightest pseudo-scalar Goldstone mass of 220 
and 260 MeV respectively [^. The calculations with the stout-link action have been performed at 
the physical value of the light quark mass. 

The purpose of this work is to investigate the quark mass dependence of the EoS by calculating 
it with the p4-action for physical values of the (degenerate) light quark masses. The calculational 
procedure used in this work closely follows that used in our previous calculations at rhi = 0.1m.s 

2. Calculations of the Equation of State 

We have performed calculations with the p4-action for fourteen values of the gauge coupling 
j3 = in the region of the finite temperature crossover. The finite temperature calculations 
have been performed on 32^ x 8 lattices, while the corresponding zero temperature calculations 
have been performed on 32^ lattices. We used the physical value for the strange quark mass and 
the ratio of strange to light quark mass was chosen to be = m^/m; = 20. The lattice spacing 
was determined by calculating the static potential and extracting the Sommer scale ro from it. 
To remove the additive divergent constant in the potential following Ref. [^] we normalized it 
to the string potential Vstnng{r) = —7i/{l2r) + ar at distance r = 1.5ro. This is needed for the 
renormalization of the Polyakov loop as discussed in Ref. [^. In Figure |l| the static potential in 
units of ro and normalized to the string potential is shown. No discretization errors are visible 
in the potential. We extracted pseudo-scalar meson masses using wall sources in the calculation 
of meson propagators. It turned out that the rjss mass is, with 1 -2 % accuracy, the same as in 
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Figure 1: The static potential in units of calculated at different gauge couplings. The dashed line shows 
the string potential. 



[^]. Thus, a re-adjustment of the line of constant physical Tj^^- mass corresponding to our new 
and smaller light quark masses was not necessary.In fact, in the present calculations we find that 
our quark mass values define a line of constant physics characterized by the following relations 
''0 • wi;r = 0.371(3), ro • niK = 1.158(5), tq • triri^- = 1.578(7) Using ro = 0.469 fm, as determined in 
Ref.[[l^], we get nijt = 154 MeV, mk = 486 MeV and^ m^,^- = 663 MeV This means that both the 
light quark masses and the strange quark mass are very close to their physical values. Furthermore, 
in the entire parameter range covered by our thermodynamic calculations deviations of the meson 
masses from the above values are less than 3%. 

The calculation of the EoS starts with the evaluation of the trace anomaly, i.e. the trace of 
the energy-momentum tensor &^^{T). It is related to the temperature derivative of the pressure 
through thermodynamic identities, 

®nniT) ^ £-3j^ d / p\ 

\ ) ' \ ■ ) 

The trace anomaly can be expressed in terms of the expectation values of quark condensates and 
the gluon action density, see e.g. Ref. [0]. The numerical results are shown in Figure ^ and are 
compared to the previous calculation at twice larger quark mass mi = 0. lm.v on A^^- = 6 lattices 
and N-^ = S lattices The differences between N-^ = 6 and Nt = S calculations are due to cutoff 
effects and have been discussed in Ref. [0].As one can see from the figure the main differences to 
the A^t- = 8 results at mi = O.lmg arise for temperatures T < 200 MeV. These differences can be 
understood as resulting from an expected shift of the transition temperature by 5 MeV when the 
light quark mass is lowered to approximately its physical value. At lower temperatures it also is 
expected that the trace anomaly increases with decreasing quark masses as hadrons become lighter 
when the quark mass is decreased. While a tendency for such an increase may be indicated by the 
data at the lowest two temperatures reached in our calculation, this effect is certainly not significant 
within the current statistical accuracy. 



A physical value for the rjss mass can be obtained from the relation ra,)^. = i / 2mj^ — — 686 MeV. 
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Figure 2: The trace anomaly {e — 3p)/T^ calculated for the physical quark mass and compared with previ- 
ous calculations at larger light quark masses m/ = Q.lnis (left) and the pressure as well as the energy density 
(right). We compare the trace anomaly with as well as with the HRG model, which includes all the reso- 
nances up to 2.5GeV. The horizontal band shows the expected uncertainty in the energy density due to the 
choice of the lower integration limit 



At temperatures below the transition temperature it is expected that thermodynamic quantities 
are well described by a hadron resonance gas (HRG) model. In fact, the freeze-out of hadrons in 
heavy ion experiments takes place in the transition region and the observed particle abundances are 
well described by the HRG model [12, 13]. Therefore in Figure ^ we also show the prediction of the 
HRG model, which includes all the known resonances up to the mass M^ax = 2.5 GeV. The lattice 
data for £ — 3/7 are below the HRG prediction although the deviations from it are smaller compared 
to the results obtained at m/ = O.lmv. We mention again the present statistical accuracy and the 
possibility of discretization effects in the hadron spectrum. In particular, due to taste breaking of 
staggered fermions pseudo-scalar mesons are not degenerate at finite lattice spacing, therefore their 
contribution to thermodynamic quantities maybe suppressed. 

From the trace anomaly the pressure and thus other thermodynamic quantities can be calcu- 
lated by performing the integration over the temperature 

^-^ = //^'^e„„(r'). (2.2, 

Here Tq is an arbitrary temperature value that is usually chosen in the low temperature regime where 
the pressure and other thermodynamical quantities are suppressed exponentially by Boltzmann 
factors associated with the lightest hadronic states, i.e. the pions. Energy e and entropy {sT = 
(p + e)) densities are then obtained by combining results for p/T^ and {e — 3p)/T'^. The numerical 
results for the pressure and energy density are shown in Fig. |2[ The uncertainties from the choice of 
the lower integration limit are shown as a horizontal band in the figure. The estimated uncertainties 
are about 8% in the energy density at the highest temperature of T ~ 260MeV, and about 13% for 
the pressure. 



3. Deconfinement and chiral aspects of the QCD transition 

In the previous section we have seen that the energy density shows a rapid rise in the tem- 
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perature interval T = (170 — 200)MeV. This is usually interpreted to be due to deconfinement, i.e. 
liberation of many new degrees of freedom. For sufficiently large quark mass this transition is 



known to be a first order transition (see e.g. Ref. [|16|]). In the limit of infinitely large quark mass 
the order parameter for the deconfinement phase transition is the Polyakov loop. After renormaliza- 
tion it can be related to the free energy of a static quark anti-quark pair Foc{T) at infinite separation 



[17, 18] Lren{T) = exp{—Foc{T)/{2T)) A rapid change in this quantity is indicative for deconfine- 
ment also in the presence of light quarks. In the opposite limit of zero quark mass one expects 
a chiral transition and the corresponding order parameter is the chiral condensate. For a genuine 
phase transition, i.e. in the chiral limit the chiral condensate vanishes at the critical temperature 
Tc. However, we expect that even for the crossover at finite quark mass the light quark condensate 
rapidly drops in the transition region, indicating an approximate restoration of the chiral symmetry. 
At non-vanishing quark mass the chiral condensate needs additive and multiplicative renormaliza- 
tion. Therefore, following Ref. 0] we introduce the so-called subtracted chiral condensate 

Here the subscripts / and s refer to light and strange chiral condensates, while the subscript and 
T to the case of zero and finite temperature respectively. Subtraction of the strange quark conden- 
sate multiplied by the ratio of the light to strange quark mass removes the quadratic divergence 
proportional to the quark mass. 

In Figure |3| we show the renormalized Polyakov loop and the subtracted chiral condensate A/ 
and compare with previous calculations performed at light quark masses equal to one tenth of the 
strange quark mass [^. The renormalized Polyakov loop rises in the temperature interval T = 
(170 — 200) MeV where we also see the rapid increase of the energy density. At the same time the 
subtracted chiral condensate rapidly drops in the transition region, indicating that the approximate 
restoration of the chiral symmetry happens in the same temperature interval as deconfinement. 
Compared to the calculation performed at light quark masses equal to one tenth of the strange 
quark mass we see a shift of the transition region by roughly 5 MeV. We note that such a shift arises 
differently in different observables. In the case of the subtracted chiral condensate, for instance, a 
major ingredient to the 'shift' is the fact, that at fixed temperature the condensate in the transition 



region is strongly quark mass dependent and drops proportional to y^mi/nis 1 19]. 

The fluctuation of strangeness is also indicative of deconfinement. It can be defined as the 
second derivative of the free energy density with respect to the strange quark chemical potential 

Z.(r) = ^ Im,^o. (3.2) 

At low temperatures strangeness is carried by massive hadrons and therefore strangeness fluctua- 
tions are suppressed. At high temperatures strangeness is carried by quarks and the effect of the 
strange quark mass is small. Therefore strangeness fluctuations are not suppressed at high tem- 
peratures. As discussed in Ref. [^] strangeness fluctuations behave like the energy density in the 
transition region, i.e. they rapidly rise in a narrow temperature interval. In Fig. ^ we show the 
strangeness fluctuations calculated at m/ = O.OSm^ and compare them with previous calculations 
performed at mi = O.lm^ M\. In the bottom figure we also show the strangeness fluctuation for 
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Figure 3: The renormalized Polyakov loop (left) and the subtracted chiral condensate (right) as function of 
the temperature calculated at m/ ~ O.OSm, and at 0.1 w^-. 
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Figure 4: Strangeness fluctuations as function of the temperature calculated at m/ = COSm^ and at O.lmj. 
In the right plot the numerical data for m/ = 0. Imj have been shifted by 5 MeV. 



mi = O.lm.5 with a 5 MeV shift of the temperature scale. As one can see this shift accounts for 
most of the quark mass dependence of the strangeness fluctuations. This is consistent with the 
conclusion obtained from the quark mass dependence of other thermodynamic observables. 

4. Conclusion 

We have calculated the EoS, renormalized Polyakov loop, subtracted chiral condensate and 
strangeness fluctuations in (2+l)-flavor QCD in the crossover region from low to high tempera- 
tures using the improved p4 staggered fermion formulation on lattices with temporal extent N-c = S 
at physical values of the light and strange quark masses. We found that thermodynamic quantities 
below the deconfinement transition are larger compared to the previous calculations performed at 
twice larger quark mass but fall below the resonance gas model result. The differences in the ther- 



6 



On Equation of State 



Peter Petreczky (for RBC-Bielefeld Collaboration) 



modynamic quantities calculated at mi = 0.05ms and mi = O-lm^ can be well understood in terms 
of the shift of the transition temperatures towards smaller values when the quark mass is decreased. 
This conclusion is also supported by the calculation of renormalized Polyakov loop, subtracted 
chiral condensate and strangeness fluctuations. No additional enhancement of the pressure and the 
energy density is seen at low temperatures. This and the deviation from the resonance gas model 
may be a cutoff effect due to taste violations. However, better statistical accuracy and calculations 
at smaller lattice spacing are needed to quantify this assertion. The transition region in our calcula- 
tions corresponds to larger temperatures compared to recent calculations with stout action |^, 
It remains to be seen whether the taste symmetry violations which are larger for the p4 action are 
responsible for this discrepancy. At temperatures above 200 MeV no quark mass dependence is 
seen in the equation of state. 
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